arXiv:1501.03798v2 [astro-ph.EP] 24 Feb 2015 


Draft version February 25, 2015 

Preprint typeset using DTeX style emulateapj v. 5/2/11 


K2-3: A NEARBY M STAR WITH THREE TRANSITING SUPER-EARTHS DISCOVERED BY K2 

Ian J. M. Crossfield^’^^, Erik Petigura^, Joshua E. Schlieder^’^^, Andrew W. Howard^, B.J. Fulton^, Kimberly 
M. Aller^, David R. CiardP, Sebastien Lepine®, Thomas Barclay^, Imke de Pater^, Katherine de Kleer^, Elisa 
V. Quintana^, Jessie L. Christiansen^, Eddie Schlafly^, Lisa Kaltenegger^\ Justin R. Crepp®, Thomas 
Henning^, Christian Obermeier^, Niall Deagon^, Lauren M. Weiss^, Howard T. Isaagson^, Brad M. S. Hansen^°, 
Mighael C. Liu^, Tom Greene^, Steve B. Howell^, Travis Barman\ Christoph MordasinT 

Draft version February 25, 2015 

ABSTRACT 

Small, cool planets represent the typical end-products of planetary formation. Studying the archi¬ 
tectures of these systems, measuring planet masses and radii, and observing these planets’ atmospheres 
during transit directly informs theories of planet assembly, migration, and evolution. Here we report 
the discovery of three small planets orbiting a bright {Ks = 8.6 mag) MO dwarf using data collected as 
part of K2, the new ecliptic survey using the re-purposed Kepler spacecraft. Stellar spectroscopy and 
K2 photometry indicate that the system hosts three transiting planets with radii 1.5 - 2.1 R 0 , strad¬ 
dling the transition region between rocky and increasingly volatile-dominated compositions. With 
orbital periods of 10-45 days the planets receive just 1.5-10 x the flux incident on Earth, making 
these some of the coolest small planets known orbiting a nearby star; planet d is located near the 
inner edge of the system’s habitable zone. The bright, low-mass star makes this system an excellent 
laboratory to determine the planets’ masses via Doppler spectroscopy and to constrain their atmo¬ 
spheric compositions via transit spectroscopy. This discovery demonstrates the ability of K2 and 
future space-based transit searches to find many fascinating objects of interest. 

Subjeet headings: K2-3— techniques: photometric — techniques: spectroscopic — eclipses 


1. INTRODUCTION 


Surveys for new planets demonstrate that sm all, low- 


mass planets are co mmon around EGK s tars (Howard 
et al.| 2QlQ, 2012). Petigura et al. (2013) usea Kepler 


data to measure the frequency of Earth-sized planets in 
Earth-like orbits to be 5-20%. Such small planets with 
moderate insolation levels (the stellar energy received by 
the planet at the top of any atmosphere) are of consid¬ 
erable interest for their ability to host Earth-like atmo¬ 
spheres that could potentially support life. 

M dwarfs offer a shortcut to observing rocky and po¬ 
tentially habitable planets. Gompared to nearby Sunlike 
stars, planets around M dwarfs are easier to find with 
transits or radial velocities (RV), they occur more fre- 
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quently (Howard et al.||2012), an d their atmospheres are 
easier to study when transiting (Stevenson et al. 2010 
Kreidberg et al.|2014). Planets transiting IVI dwarfs otter 
the best opportunity to study habitability and con strain 
models of rocky planet a ssembly and migration (Swift 
et al.||2013| |Hansen||2014|) and of planetary atmosphere s 
et al. |20Il|[Rodler fc LQpez-Morales|2014[) 
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ulti-planet IVI dwarf systems are even more exciting, 

both because such candidates are extre mely unlikely to 
result from astrophysical false positives (Lissauer et al. 
2012) and becau se they allow for studi es of compara- 
tive planetology (Muirhead et al. 2012) with identical 
initial conditions (i.e., formation in the same natal disk). 
However, relatively few confirmed transiting planets (and 
fewer multiple systems) are known around M dwarfs, 
and the (because Kep/er’s prime mission targeted just 
3900 late-type dwarfs) the prevalence of planets around 
M dw arfs is less well constrained th an around Sunlike 
stars (Dressing & Gharbonneau||2013). 

We are using K 2, the contihiimg m ission of NASA’s 
Kepler spacecraft (Howell et al.|2014), to target thou¬ 
sands of M dwarfs in each K2 held to find new, small 
planets orbiting these stars. K2’s 80-day campaigns 
are ideally suited to finding large numbers of small, 
cool planets around M dwarfs, out to semimajor axes 
in th e stars’ habitable zones (e.g., Kopparapu et al. 


2014). In addition, some of K2’s M-dwart planets orbit 


stars bright enough for atmospheric character ization via 
JWST transmission or emission spectroscopy (Kalteneg- 
ger & Traub||2009t |Batalha et al.||2013| IBeichman et al 

Here, we present the discovery of a new multi-planet 
system orbiting a bright M dwarf (K2-3, PMIl 1293-0127, 
UGAG4 443-054906, PPMX 112920.3-012717). We de¬ 
scribe our analysis of the K2 photometry and of supple- 
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TABLE 1 

Stellar Parameters of K2-3 


Parameter 


Value 


Source 


Identifying information 


a R.A. (hh:mm:ss) 

11:29:20.388 



(5 Dec. (dd:mm:ss) 

-01:27:17.23 



EPIG ID 

201367065 



2MASS ID 

11292037-0127173 

2MASS 


Photometric Properties 



B (mag). 

13.52 ± 0.06 

APASS 


V (mag). 

12.17 ± 0.01 

APASS 


g (mag) . 

12.871 ± 0.030 

APASS 


r (mag) . 

11.582 ± 0.020 

APASS 


i (mag). 

10.98 ± 0.17 

APASS 


J (mag). 

9.421 ± 0.027 

2MASS 


H (mag) . 

8.805 ± 0.044 

2MASS 


Ks (mag) . 

8.561 ± 0.023 

2MASS 


W1 (mag) . 

8.443 ± 0.022 

AllWISE 


W2 (mag) . 

8.424 ± 0.019 

AllWISE 


W3 (mag) . 

8.322 ± 0.021 

AllWISE 


Spectroscopic and Derived Properties 



flex (mas yr 

88.3 ± 2.0 

Zacharias et al. 

(2012 

fi5 (mas yr“^) 

-73.6 ± 2.7 

Zacharias et al. 

(2012 

Barycentric rv (km s“^) 

32.6 ± 1 

APh', this paper 

Distance (pc) 

45 ±3 

this paper 


EW (Ha) (A) 

0.38 ± 0.06 

EEOSG, this paper 

Age (Gyr) 


EFOSG, this paper 

Spectral Type 

M0.0T0.5V 

This paper. 


[Fe/H] 

-0.32 ± 0.13 

SpeX, this paper 

Teff (K) 

3896 ± 189 

SpeX, this paper 

M* (Mq) 

0.601 ± 0.089 

SpeX, this paper 

R* (Rq) 

0.561 ± 0.068 

SpeX, this paper 


merit ary imaging and spectroscopic data in Sec. i In 
Sec. we present the results of our analysis of K2-3’s 
properties and discuss the potential for future observa¬ 
tions of this and other systems discovered by K2. 


2. OBSERVATIONS AND ANALYSIS 

We identified the high proper motion star PMI11293- 
0127 as a target for our Campaign 1 proposal (G01Q3q^ 
PI Cros sfield) from the SUPERBLINK proper motion 


survey (Lepine & Shara 2005 Lepine fc Gaidos| 


losl|mT|. 

"M dwarf 


We identified the star as a probable nearby 
based on a color and proper motion selection scheme and 
selecting all targets with {V — J) > 2.5, ^ + 5 log /i + 5 < 
10, and {6V — 7J — 3) < 51og/i, where /i is the proper 
motion. The s tar matched t he source K2-3 in the Kepler 
input catalog ( |Huber|[2014 ). K2 then observed this tar¬ 
get in long-cadence mode during Gl, covering 30 May 
to 21 Aug 2014. Target properti es, including optica l 
and NIR photometry from AP ASS dHende n et al.||2012|), 
2MA SS dSkrutskie et al.|2006|, and WT5E7[Wngnt^r^ 
2010) are summarized in 'Table [^. 


2.1. K2 Photometry 

2 . 1 . 1 . Extracting the Photometry 

We extracted the photometry K2-3 from the pixel data, 
which we downloaded from the MAST. Because K2 only 
has two functional reaction wheels, the telescope cannot 
maintain the 50-milhpixel pointing precision achieved 
during the prime mission. The dominant drift is roll 

The star was also identified in programs GO1006, GOI050, 
GO1052, GO1036, GOI075, GOI059, and GOI063. 


around the telescope boresight. When the spacecraft 
reaches a pre-determined limit the spacecraft corrects 
this roll with a thruster fire. As the spacecraft rolls, stars 
move over different pixels having different sensitivities. 
Thus, motion of the star results in apparent changes in 
stellar brightness. 

Because a target star traces out similar paths during 
each roll of the spacecraft, it is possible to separate out 
variations in stellar brightness that are roll angle depen¬ 
dent, and to remove these variations from the photome- 
try. Our extraction pipe l ine dr aws heavily on the work of 
Vanderburg & Johnson (2014). We begin by computing 
the median flux tor each trarne and adopt this value as 
the background flux level. The background flux is sub¬ 
tracted out on a frame by frame basis. We compute the 
raw photometry, Fsap, by summing the flux within a 
soft-edged circular aperture centered around the target 
star. We compute the row and column centroids within 
the aperture. 

On short timescales, spacecraft roll is the dominant 
motion term and can be described by a single variable. 
We identify the roll direction by computing the principle 
components of the row and column centroids, x' and y'. 
We fit for a function that relates Fsap to x'. We describe 
this trend by Fsap = GP(t'), where GP is a Gaussian 
process having a correlation matrix given by a squared 
exponential kernel. Fitting the GP(t') is an iterative 
process where outliers are identified and removed and 
the hyperparameters associated with the squared expo¬ 
nential kernel are adjusted to yield the minimum residual 
RMS. 


The algorithm described in Vanderburg fc Johnson 


(2014) was developed for the K2 engineering campaign 
(CO), where the time baseline was short enough that 
drifts in stellar position along the y' direction could be 
ignored. During 80 day period of Gl observations, stars 
moved enough along the y' direction that the GP(t') de¬ 
termined using data early in the campaign was no-longer 
an appropriate description of the position-depen dent flux 
varia t ions. Adopting an approach described in [Vander^ 
burg| (|2014|), we divided the Gl observations in to six 
nearly equal segments and performed the ID decorrela¬ 
tion approach described above on each segment individ¬ 
ually. The entire procedure described above is repeated 
for different aperture radii (2, 3, 4, 5, 6, and 7 pixels). 
We select the aperture size that minimizes the calibrated 
RMS. For K2-3, a circular aperture with a 4 pixel radius 
yielded the best calibrated photometry (which is avail¬ 
able as an electronic supplement to this paper). Our 
photometry suggests that K2-3 may exhibit photomet¬ 
ric variations of < 1% on week-to-month timescales, but 
K2’s long-term stability is constrained sufficiently poorly 
that we cannot claim evidence for periodic modulation 
indicative stellar rotation. 

2 . 1 . 2 . Transit Detection 

We searched through the calibrated and detrended 
photometry (show n in Fig. us ing th e TERRA algo¬ 
rithm described in|Petigura et al.| (|2C)13|). TERRA iden¬ 
tified a transit candidate having E = 10.056 days and 
SNR = 59. We fit this candidate with a [Mandel fc Agol| 
(2002) model and subtracted the best fit model from the 
photometry. We reran TERRA on the photometry with 
the P = 10.056 day candidate removed. We found a sec- 
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ond candidate having P = 24.641 days and SNR = 30. 
Again we removed the best-fitting model. TERRA did 
not find any additional transits, but a ^45-day candidate 
was identified by eye (TERRA currently requires 3 de¬ 
tected transits, and thus was not sensitive to the longest 
period candidate which only transits twice during Cl). 
We fit each of these two transits individually and find 
consistent transit parameters, supporting the hypothesis 
that they result from a single planet. At half of this pe¬ 
riod a third transit would occur in Cl’s data gap (see 
Fig. 0, but this would give the outer two planets a pe¬ 
riod ratio of just 1.1. The previous recor d-holder for a 
close period ratio is the Kepl er-36 system (Carter et al. 
2012| I Winn & Eabrycky|2014|) , whose two planets exhibit 
a considerably larger period ratio of 1.17 and transit tim¬ 
ing variations of many hours. It is unlikely that such an 
unusual system would he just 45 pc away; in addition, our 
dynamical analysis (described below) indicates that this 
period ratio would be dynamically unstable. We there¬ 
fore conclude that the third planet’s period is ^45 d. 

2.2. Target Validation and Stellar Speetroseopy 

We conducted a number of pixel-level diagnostics and 
observed K2-3 using several spectrographs to constrain 
the stellar properties. These observations are described 
below. The reduced spectra are attached as an electronic 
supplement to this paper, and the derived parameters are 
listed in Table [H 

2 . 2 . 1 . Pixel-Level and Photometrie Data Validation 

Experience over the last decade shows that transit-like 
signals must be validated to ensure that they arise from 
true planets, not “false-positive” configurations such as 
backgroun d eclipsing binaries blend ed with foreground 
stars (e.g., |Torres et al. 2004 2011). We therefore im¬ 
plement a large number of tests on the pixel-level data 
and extracted photometry to identify and weed out these 
false positives. 

Once transit-like events are identified, TERRA runs a 
suite of diagnostics to distinguish planets from phenom¬ 
ena like eclipsing binaries, starspots, and other periodic 
stellar variability. We subject targets passing this first 
step to an extensive battery of further tests which search 
for blends using an examination of centroid motions in 
and out of transit, difference imagin g analyses, and con - 
struction of pixel correlation images ( Bryson et al.|2Q13 ). 
Though we are still learning how to optimally tune these 
tests to account for K2’s few-pixel pointing variations, 
validation results for large numbers of targets indicates 
that the transit-like events identified with K2-3 occur 
within roughly one pixel (4”) of the target star. When 
combined with our seeing-limited and ada ptive optics 
imaging described below, as described in Sec. |2.5| we find 
that K2-3’s transits are far more likely to be explained by 
a multi-planet system than by nonplanetary phenomena. 

2 . 2 . 2 . Optieal and Infrared Speetroseopy 

We obtained R ^ 1500 spectra from 0.6-1.0jam of K2-3 
and a number of calib ration objects using NTT/EEOSC2 
dBuzzoni et al.||1984|) on UT 11 Jan 2015 as part of 70- 
night K2 followup program (BID 194.C-0443, PI Cross- 
field). We draw our ca l ibrators from several recent works 
(Boyajian et al. 2012 Recant & Mamajek 2013 Mann 


et al.||2013a). A forthcoming paper will discuss these ef¬ 
forts; uTTnef, we bias-subtract and fiat-field the data 
frames, extract spectra using IRAE, and wavelength- 
calibrate using EEOSC2’s internal HeAr lamps. We 
achieve a S/N per resolution element of ^100 for K2- 
3 and somewhat higher for our reference sample. We 
flux-calibrate the extracted spectrum using observations 
of spectrophotometric standards. 

We observed K2-3 on 2015 January 11 UT using the 
recen tly refurbished SpeX spectrograph (Rayner et al. 
20031) on the 3.0m NASA Infrared Telescope facility 
(IRl E). The data were taken under clear skies with an 
average seeing of ^077. We observed with the instrument 
in short cross dispersed mode (SXD) using the 0.3 X 15” 
slit. This setup provides simultaneous wavelength cover¬ 
age from 0.7 to 2.5 jam at a resolution of R?^2000. The 
extended blue wavelength coverage is a result of the re¬ 
cent chip upgrade SpeX received in July 2014. The target 
was placed at two positions along the slit and observed in 
an ABBA pattern for subsequent sky subtraction. The 
observing sequence consisted of 8 x 40 s exposures for a 
total integration time of 320s. Once the exposures were 
stacked, this integration time led to a signal-to-noise of 
> 140 per resolution element. We obtained standard 
SpeX calibration frames consisting of fiats and arclamp 
exposures immediately before observing K2-3. 

The SpeX spectrum was reduced using the SpeX- 
Tool software package (jCushing et al.||2004). SpeXTool 
performs fiat-field correction and wavelengtn calibration 
from the calibration frames followed by sky subtraction 
and extraction of the one-dimensional spectrum. Indi¬ 
vidual exposures of the target were combined using the 
xeomhspee routine within SpeXTool. We corrected for 
atmospheric absorption and performed flux calibration 
using the AOV-type star HD 97585 which was observed 
within 20 minutes and 0.015 airmass of the target. A 
telluric correction spectrum was constructed fro m the 
spectrum o f the AOV using the xtelleor package (Vacca 
et al. 2003) and applied to the spectrum of K2-3. This 
package also performs flux calibration. Separate, telluric- 
corrected SpeX orders were combined and flux matched 
into a continuous spectrum using the xmergeorders rou¬ 
tine. To minimize errors in the spectral slope due to 
changes in seeing, guiding, and differential refraction, we 
aligned the slit with the parallactic angle and minimized 
the time between observations of the target and standard 
star. Prior to performing any spectroscopic analyses, we 
also applied corrections for the barycentric velocity of 
the observatory and the measured radial velocity. The 
final, calibrated spectrum is shown in Eig. 

2 . 2 . 3 . Stellar Parameters 


Mann et al.|( 2013b) motivate a set of temperature sen- 
sitive spectral indices spanning the visible, J-, H-, and 
K-bands that a r e calib rated using the M dwarf sample of 
Boyajian et al. (2012) with interferometrically measured 
radii. We used these indices to estimate the temperature 
of K2-3. We calculate the mean of the temperatures from 
each of the three NIR band indices and their rms scatter 
and find Teff = 3896T117 K (±148 K systematic error, 
±189 K total error). This range of effective temperatures 
is consistent with that main-sequence MO dwarfs o f spec¬ 
tral type K8V to MOV rtPecaut & Mamajek|2013|) , using 
the modified system which incorporates subtypes K8V 
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Fig. 1.— Top: Calibrated K2 photometry for K2-3. Vertical ticks indicate the locations of each planets’ transits. Bottom: Phase-folded 
photometry and best-fit light curves for each planet. 
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Fig. 2. — Calibrated IRTF/SpeX spectra of our target compared to spectral standards. Stellar parameters are tabulated in Table 


and K9V between K7V and MOV. 

We a dopt the metallicity calibration of [Mann et al. 


(2013a) to remain consistent with our methods for de- 


termmmg Teff, and other parameters. We use custom 
IDL software provided by A. Mamf^ to calculate the 
metallicity in in the vi sible, J-, H-, and K-b ands follow¬ 


ing the calibrations of Mann et al. (2013a). Since our 


SpeX spectrum does not extend <0.7 /im, we do not use 
the visible band cali brations. Following the suggestion of 
Mann et al. (2013a), we also discard the J-band metal- 
licity, which is often an outlier. Our final metallicity 
is the mean of those measured from the H- and K-band 
relations and the error is the quadrature sum of the mea¬ 
surement error and systematic error in each band. We 
find [Fe/H] = -0.32±0.13. Thus, K2-3’s metallicity is 
sub-solar, broadly consistent with many other nearby, 
fi eld-age, M dwarfs. 


Mann et al. (2013b) provide empirical calibrations to 
calculate the radii, niasses, and luminosities given the 


Teff of an M dwarf. We estimate these additional fun¬ 
damental parameters again using IDL software written 
by A. Manrp^ to calculate radius, mass, and luminosity 
and their associated errors using the relations detailed in 
Using the most conservative Teff 
Rq 


Mann et al. (2013a). 

culate = 


errors, we calculate = 0.561 ±0.068 Rq and = 
0.601±0.089 Mq. These values, and the other funda¬ 
mental parameters of the star, are tabulated in Table 
and are used for subsequent estimates of the individual 
planet properties. 

Independent of these parameters, we also assign a spec¬ 
tral type to this star using molecular band heads in our 
optical and NIR sp ectra. In the optical, the Ti0 5, CaH2, 
and CaH3 indices ( |Reid et ^|1995| [Giz^[j ^"7|) are cali¬ 
brated for the earliest IVl dwarfs ( [Lepine et al.|2003 ) and 
avoid regions of the spectrum with heavy telluric contam- 
ination. Following the most recent spectral type calibra¬ 
tions of these indices by Lepine et al. (2013|), our EFOSC 
spectrum yields a spectral type of K7.5±0.5, determined 


https://github.com/awmann/metal 


https://github.com/awmann/Teff_rad_mass_lum 
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to a half-subtype scale and assuming a sequence K5-K7- 


MQ (i.e., without the K8 and K9 subdivisions of Pecaut 
& Mamajek||2013j). In the NIR, the H 2 O-K 2 index mea¬ 
sures water opacity in t he K-band, and wa s calib rated 
to a spectral subtype by Rojas-Ayala et ah (2012). We 
calculate this index fromi our bpe X spectriirn and esti- 
mate a spectral type of M0.5±0.5. Lepine et ah (2013) 
also provide a calibration of the V — J color to spectral 
subtype. Our target has V — J = 2.75, consistent with 
subtype K7.5 on the scale of Pecaut & Mamajek (2013). 

The spectroscopic and photometric classitications 
are all consistent, although the NIR classification is 
marginally later. Here we average the optical and in¬ 
frared results and adopt a spectral type of M 0.0T0.5V. 
Using the riJi7A" photometric calibrations of |Krans~S^ 
Hillenbrand (2007), we estimate a distance to K2-3 of 
45 ± 3 pc. - 

We obtained high-resolution (2” slit width with the 
B decker) spectra of K2-3 with the Levy Spectrometer 
([Radovan et al.|| 2010|) on the Aut omated Planet Finder 
(aPI^ 4 telescope ([Vogt et al.|2014 ). The spectra were re- 


dnced usin g standard procedures, as described in (Fulton 
et al.|20f5 ). Inspection of the gravity-sensitive lines 


1 con- 


hrms that K2-3 is a high gravity target, consistent with 
the medium resolution spectra described above. We do 
not see any evidence of a second set of spectral lines, 
ruling out companions ^2.5 mag fainter than K2-3 at 
visible wavelengths. 

2.2.4. Activity, Age, and Membership 

Lines in the Balmer series are associated with mag¬ 
netic activity in late-type stars. The strongest line in 
the series, Ha at 6563 A, is classically used to ass es the 
activity of M dwa rfs and as a crude indicator of age (West 
1|2QQ4| 

] width 


et al. 


alent 


2008). We therefore measure the Ha e qmv- 
(FW) as defined by West et al. (2011) and 


Lepine et al. (2013) and find consistent results using both 
approaches. We use two different integration regions to 
calculate this EW and apply Monte Carlo methods to 
estimate the uncertainty in the EW measurements. We 
find that EW = 0.38 ± O.OOA, indicating that K2-3 is a 
relatively inactive star. We further investigate possible 
chromospheric activity in K2-3 by analyzing its UV emis¬ 
sion measured by GALEX ([Martin et al.|2005|). The star 
is a weak near-UV (NUV) emitter and is not detected in 
the far-UV (FUV). Its low NUV flux and non-detection 
in the FUV is consistent with qui escent emission, simi - 
lar to other nearby field M dwarfs (Shkolnik et al.|2011). 
The Ha absorption, UV fluxes, and lack of chromospheric 
activity in an MO dwarf all indicate an old, fi eld-age star 
and t ranslates to a lower age limit of ^1 Gyr (West et al. 
2008|. 

We further examined the possibility that K2-3 is young 
by comparing its space position {XYZ) and kinematics 
(UVW) with those of known young moving groups. Its 6- 
dimensional UVWXYZ position is inconsistent with the 
well-kno wn nearby young (^ 10-100 Myr) groups summa¬ 
rized by Gagne et al. ( 2014[), as well as other sparser or 


et al.|2013[). To provide a quantitative estiniate, we used 


slightly old er groups (bhkolnik et al. [2009 Zuckerman 


the B ANyAN II web tool (Malo et al.|2013l|Gagne et al. 
2014). BANYAN II calculates the probability of an ob¬ 
ject being a member of a nearby young moving group 


using Bayesian inference and the observed proper mo¬ 
tion, sky coordinates, radial velocity and distance. The 
probability of K2-3 being a member of one of the known 
nearby young moving groups is <0.1% given the sky co¬ 
ordinates, proper motion, and radial velocity. Inclusion 
of the photometric distance estimate (and conservatively 
assuming a 20% distance uncertainty) does not change 
the BANYAN II results. Thus we conclude that K2-3 is 
unlikely to be a member of any of these young moving 
groups. 

2.3. Archival and Adaptive Optics Imaging 

To rule out the presence of a background star being the 
source of or diluting the transit events, we compare two 
epochs of imaging data from the Digitized Sky S urvey 


(DSS) and the Sloan Digital Sky Survey (SDSS; [Ahn 
et al.[[2012[) separated by 45 years. The data shown 


in I^'ig. are the DSS-Red plates with a pixel scale of 
1.7"/pixel taken on 19 April 1955 and the SDSS r-band 
image with a pixel scale of 0.396'Ypixel taken on 03 
March 2000. The images are 1 arcminute on a side and 
clearly show the proper motion of the primary target. 
The nearby star located 27" to the NE is consistent with 
zero motion within our astrometric uncertainties; this 
star lies outside the photometric aperture applied to the 
K2 photometry. The primary target, in contrast, dis¬ 
plays a clear proper motion of 6.2" over 45 years, in 
reasonable agreement with the measured prope r motion 
(Lepine & Gaidos 2011 Zacharias et al. 2012). In the 
DSS image there is no evidence of a background star, 
and we estimate if a star is located at the position of the 
primary target in the Kepler data, that star must be at 
least 6 magnitudes (or more) fainter than the target star. 

Near-infrared adaptive optics imaging of K2-3 was ob¬ 
tained at Keck Observatory on the nights of 2015 Jan¬ 
uary 12 UT and 2015 Janary 16 UT. Observations were 
obtained with the 1024 x 1024 NIRC2 array and the nat¬ 
ural guide star system; the target star was bright enough 
to be used as the guide star. The data were acquired in 
the narrow-band K-band and J-band continuum filters 
(Kcont and Jcont) using the narrow camera field of view 
with a pixel scale of 9.942 mas pix“^; the atmosphere 
was less stable on night 1 and only Kcont was acquired 
on that night. A 3-point dither pattern was utilized to 
avoid the noisier lower left quadrant of the NIRG2 ar¬ 
ray. For both nights, the 3-point dither pattern was ob¬ 
served with 10 coadds and a 1.5 integration time, but 
on night 1 only 4 frames were acquired for a total of 60 
seconds of on-source exposure time. For night 2, three 
full dither patterns were acquired for a total on-source 
exposure time of 135 seconds in both Kcont and Jcont 
filters. The data from each night were flatfielded and 
sky subtracted and the dither positions were shifted and 
coadded into a single final image. The final images from 
night 2 are shown in Fig. and d. 

For night 1, the target star was measured with a reso¬ 
lution of 0.07" (FWHM), but the atmosphere was much 
more stable during night 2 and these images have a reso¬ 
lution of 0.05" in the Kcont filter and 0.04" in the Jcont 
filter. No other stars were detected within the 10" field 
of view of the camera; speckles seen in the Kcont images 
are not co-spatial with the speckles seen in the Jcont im¬ 
age, indicating that the speckles are not faint compan¬ 
ions. The night 2 data were much more sensitive than 
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the night 1 data and we report the analysis of those data 
in this work. In the Kcont filter, the data are sensitive to 
stars that have K-band brightness of Amag = 2.4 mag 
at a separation of 0.05" and Amag = 8.0 mag at a sep¬ 
aration of 0.5" from the central star; in the Jcont filter, 
the data are sensitive to stars that have J-band bright¬ 
ness of Amag = 2.0 mag at a separation of 0.05" and 
Amag = 7.5 mag at a separation of 0.5" from the cen¬ 
tral star (see Fig.j^ and d). We estimate the sensitivities 
by injecting fake sources with a signal-to-noise of 5 into 
the final combined images at distances of N*FWHM from 
the central source. The 5cr sensitivities, as a function of 
radius from the star are shown in Fig. and d. 

2.4. Light Curve Fitting 

We analyze the photometry using standard Python- 
based minimizers, t he emcee Markov Chain Monte- Carlo 
(MCMC) package (Foreman-M ackey et al. '2013|), and 
the JKTEBOP li ghtcurve code ( [bouthworth et al.||2004 


.|2011 ) using numerical integration to account 
for our ^30-min cadence. We fit each planet’s transit 
separately, after first masking out data taken during the 
other planets’ transits. 

We use the best-fit TERRA parameters to initialize the 
fits. We assumed a linear limb-darkening relation for the 
star. Because the data are insufficient to break all de- 


generacies b etween the light curve parameters (Muir head 
et al.||2Q12 ), we impose Gaussian priors in our analysis. 
I^br the lirnb-darkening parameter we assume a dis¬ 
tribution with center 0.560 and dispersion 0.044; these 
values correspond to the mean and standard deviation, 
res pectively, of all line ar limb-darkening terms tabulated 
by (Claret et al.|2012) that satisfy 3300 < Tg// < 3700 K 
and iogiQg > 4.5. Using the spectroscopic parameters 
presented below (Table [R, we also i mpose a prior on 
the stellar de nsity to constrain R^/a (Seager & Mallen- 
Ornelas||2003 ). This last point assumes that the planets' 
orbits are circular, an assumption that future RV mea¬ 
surements will test. 

We seed our 60 MCMC chains with values near the 
best-fit parameters. We assign our data points equal 
weights, such that the best-fit likelihood equals —x^/2. 
After burn-in we run the MCMC sampler: after each set 
of 2000 steps, we optimize the fits given by each chain’s 
parameters to check for better fits to the data. We re¬ 
initialize the sampler and re-scale the data weights if 
we find an improved fit, repeating until all parameters’ 
chains are well- mixed (as indicated by Gelman-Rubin 
metrics < 1.03; Gelman & Rubin 1992). As our final 
confidence intervals, we use the 15.87% and 84.13% per¬ 
centiles of each parameters’ posterior distribution. The 
final distributions are unimodal. Fig. shows the result¬ 
ing photometry and best-fit models, and Table sum¬ 
marizes the final values and uncertainties. 

2.5. Ruling Out False Positives 

Almost all candidat es in Kepler^’s multi-p lanet systems 
are bona fide planets ( Lissauer et al.||201l| , but one per¬ 
nicious source of confusion is the possibility of mistaking 
blended stars each hosting their own planets for a sin¬ 
gle multi-planet system. We therefore investigated the 
possibility that K2-3 might be a blend of multiple stars. 
First, we note that K2-3’s proper motion (listed in Ta¬ 


ble [R is large enough that optical DSS survey images re- 
vealno objects at the star’s current location (see Fig. 
and b). Our Keck/NIRC2 images also show no compan¬ 
ions at s^arations down to a fraction of an arc second 
(see Fig. and d), and so our data validation tests in¬ 
dicate that the transits must occur around K2-3 and not 
around some other nearby star. Blends involving back¬ 
ground eclipsing binaries are thus strongly disfavored. 

The most likely remaining false positive configuration 
involves a heirarchical triple system, with a later-type M 
dwarf close to K2-3 and with its own transiting planet (s) 
— but this too is extremely unlikely. An M4 dwarf 
would have AKp ^ 2.7 and so might be missed in 
our APF a nd EFOSG spectra, but the M4 would have 
AKs ~ 2.0 (jKraus & Hillenbrand||2007|) and so to avoid 
detection in our Keck/N IRC2 image it would need to lie 
at a < 2.3 AU — while still needing to host its own 2R^ 
transiting planet. The likelihood that K2-3 has such a 
low-mass stellar companion is ^0.4 and that such a com- 
panio n would lie at a projecte d separation <2.3 AU is 
^0.5 (Duchene & Kraus 2013). Eor planet b, the like- 
lihood of an M dwarf hosting such a planet is < 0.15 
(Dressing & Gharbonneau|2013); and the likelihood of it 
transiting is ^0.02. 'I'hen the likelihood of such a con¬ 
trived configuration is just ^ 6 X 10 ^ (1 in 1700), so we 
eliminate this scenario as well - in any case, such a binary 
would be quickly revealed by even crude radial velocity 
measurements. We therefore conclude that K2-3 indeed 
hosts a three-planet system. 

2.6. System Stability 

Here we investigate the dynamical stability of the 
three-planet K2-3 system. The planet masses are uncon¬ 
strained by transit photometry, so we adopt the following 
mass-radius relationship: 

• M = ^R^p, where p = (2.43 + 3.39 * {Rp/R^)) g 
cm“^ for Rp < 1.5i70 ( Weiss fc Marcy|[~2014[ ) 


• M = 2.69M^ 

1.5-R0 < Rp < 4.0-/7g 


(fe) 


0.93 


(Weiss & Marcy 


2014) for 


.M = Me(|f) 


2.06 


for Rp > 4 .OR 0 (Lissauer et al. 


2012 ) 


Adopting the above mass-radius relationship we derive 
masses of 5.3, 4.3, and 4.4 for planets b, c, and d 
respectively. We integrate the system forward in time 
with the Mercury integration package (Ghambers 1999) 
utilizing the hybrid integrator and found the system to 
be stable for the full 2 x 10^ yr simulation. 

We also evaluate analytically the system’s stability. 
The relevant length scale for dynamical interactions be¬ 
tween planets is the mutual Hill radius: 


Rh = 


Min + Afo 

3M, 


nl/3 


■ ^out 


( 1 ) 


where M and a denote mass and semi-major axis, re¬ 
spectively. The subscripts “in” and “out” correspo nd to 
the inner an outer planets respectively. Eollowing |Eab- 
rycky et al. (2012), for each pair of planets, we compute 
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Fig. 3.— We detect no objects within 25” of K2-3: in DSS (a); in SDSS (b); and with Keck/NIRC2 in J-band (c) and K-band (d). 



TABLE 2 

Planet Parameters 


Parameter 


Units 

b 

c 

d 

To 

BJDtdb - 2454833 

1980.4189lHo[l 

1979.27861^027 

1993.2232l°°°|^ 

P 

d 


10.05403l°;°[l°26 

24.6454l°;[l[;[3 

44.563ll°;[;°f5 

i 

deg 


89.28l°;« 

89.55l°j® 

89.68l°Ji 

Rp/R^ 

% 


q acQ+0.123 

9 7SfiT0.143 

'oe_Q Qg3 

2.48l°14 

Ti4 

hr 


9 c'c'q+0.047 
Z.OOO_Q Q44 

q A 90+0.106 

3-981^15 

R^ ! a 

- 


0 r)Q4Q+0-0049 

U.UC)TC)_q QQ2Q 

u.uiyc)_Q QQ^4 

0 Ql 97+0-0025 

b 

- 


0 47+0-22 
-0.23 

0 41+0-26 
*^•^^-0.25 

*^•^^-0.28 

u 

- 


U.ODU_q Q42 

U.OO/ _Q Q44 

U.ODC)_q Q42 

a 

AU 


U.U/Dy_Q QQ4Q 

U.iOyy_Q QQ73 

0 9Q76+0-0098 

U.ZU/ 0_Q QJ^Qg 

Rp 

R® 


^•-^^-0.26 

1 ^79+0.23 
‘ ^-0.22 

1 49+0-21 

J..o^_0 20 

‘8inc 

s® 


II 0+^-^ 

q q9+1.25 

O.C)Z_o 95 

1 41+0.57 

1-Ol_o.43 

Teq 

(K) 


-800 

-450 

-300 


A = (ttout — ^^in) the separation between the plan¬ 

ets measured in units of their mutual Hill radii. If two 
planets begin on circul ar orbits, they ar e stable indefi¬ 
nitely if A > 2\/3 « 3.5 


Gladman 


(1993). In the case of 


K2-3, Abe = 15.9 and Acd = 11.0. Thus, the two pairs 
of adjacent planets do not violate the criterion of Hill 
stability. 

There is no analytic stab ility criterion for system s hav¬ 


ing three or more p lanets Fabrycky et al. (2012). Fab¬ 
ry cky et al. (2012) introduce Ain+ Aout^ as a heuris- 


tic metric tor assessing the stability of three planets 
in triple or higher multiplicity systems. They adopt 
Ain T ^out > 18 as a heuristic criterion for the stability 
of three planets, n iotivated by suites of dire ct numerical 
integrations (e.g. |Smith & Lissauer (2009|)). This cri¬ 
terion is empirically supported by the ensemble of sys¬ 
tems with three or more transiting planets fro m the Ke¬ 
pler mission . Among the 413 such systems in |Fabrycky| 
et al.| (| 2012 |), only six had Ain + Aout < 18. For K2-3, 
Abe + Acd = 26.9, and thus has a similar architecture 
to the ensemble of triple and higher systems discovered 
during the prime Kepler mission. 

3. DISCUSSION 

Our analysis indicates three small planets orbiting this 
bright, nearby M dwarf. The planets range in size from 
2/^0 to 1 . 5 i? 0 , indicating that they may span the gap 
between rock-dominated “Earths” / “super-Earths” and 
low-dens ity “sub-Neptunes” with considerable volatile 
content (|Marcy et al.||2014| |Rogers||2014| jDressing et al. 
M4|. 

The planets’ radii imply masses of roughly 4-5 Me 
and Doppler amplitudes of 1.2-2.3 ms“^, within reach 


of modern RV spectrographs. These mass estimates 
assume that the planets fall on the mean mass-radius 
relationship, characterized by high densities and rocky 
compositions for planets smaller than ^1.6 Re- How¬ 
ever, most of the planets with measure d masses and 
Rr, < ^QRe have high inc i dent fluxes (e.g.j Batalha et al. 


2011 Howard et al. 2013 Pepe et al. 20131). The mass- 
radius relationship is poorly constrained tor cool planets 
that are less likely to b e sculpted by therm al evolution 
and photo-evaporation (Lopez et al. 2012 ). Character¬ 
izing the mass-radius relationship lor these cool, small 
planets is an important step to learning whether Earth- 
size planets in the habitable zone also have Earth-like 
atmospheres. 

The planets’ receive insolation levels (5'inc) roughly 11, 
3.3, and 1.5 x that of the Earth for planets b, c, and d, 
respectively. Planet d is located at the inner edge of the 
system’s habitable zone, with 5'in 


1 . 5 llo. 475 'ffi - close 


to the limits of the empirical habitable zone (e.g., Kop 
parapn et al.pofd )- making this planet a very interest- 


ing potential super-Venus or super-Earth. Because this 
system is so close the atmosphere of this planet can be 
explored in the near future; depending on atmospheric, 
cloud, and surface prope rties liquid water could poten¬ 
tially pers ist on planet c ( Zsom et al.|2Q13 but see Kast- 


ing et al. 2014) 

The K2-3 system is a convenient system to measure 
the atmospheric properties of small, cool planets. In¬ 
deed ^_the_sta£ 2 ^_^ullmagnitude brighter than Kepler- 
138 (Kipping et al. 2014), the previous best system for 
characterizing cool, nearly Earth-size planets. For cloud- 
free, hydrogen-dominated atmospheres, we estimate that 
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these planets will show spectral features wit h amplitudes 
oflOHR^jRl on the order of 100-200 ppm (Miller-Ricci 


et al. 2009), where H is the atmospheric scale height. 


These features would be detectable with current instru- 


mentation o n the Hubble Space Telescope (Kreidberg 
2014| ). Transit features in a heavy atmosphere 


et al. 


(e.g., 1^2, CO2) would be an order of magnitude smaller, 
and secondary eclipses will have depths on the order of 
50-150 ppm - either of these scenar¬ 
ios should be detectable with JWST. By allowing us to 
measure masses and atmospheric conditions for 3 small 
planets in a single system, K2-3 represents an exciting 
opportunity to test theories of planet formation and evo¬ 
lution in a single extrasolar laboratory. 

That K2 should reveal such a system in its first full 
campaign demonstrates that the mission will extend 
Kepler’s compelling scien tific legacy for years to c ome. 
Along with HIP 116454 (Vanderburg et al. 2014), the 
discovery of K2-3 shows that K2 is already finding fasci¬ 
nating new targets for observation with JWST and her¬ 
alds an era of further unprecedented discoveries in the 
TESS era. 
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